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Abstract 


Peccei-Quinn symmetry with attendant axion is a most compelling, and perhaps the 
most minimal, extension of the standard model, as it provides a very elegant solution to 
the nagging strong CP-problem associated with the 0 vacuum structure of QCD. However, 
particle physics gives little guidance as to the axion mass; a priori , the plausible values 
span the range: 10 -12 eV < m 0 < 10® eV, some 18 orders-of-magnitude. Laboratory 
experiments have excluded masses greater than 10 4 eV, leaving unprobed some 16 orders- 
of-magnitude. Axions have a host of interesting astrophysical and cosmological effects, 
including, modifying the evolution of stars of all types (our sun, red giants, white dwarfs, 
and neutron stars), contributing significantly to the mass density of the Universe today, 
and producing detectable line radiation through the decays of relic axions. Consideration 
of these effects has probed 14 orders-of-magnitude in axion mass, and has left open only 
two windows for further exploration: 10“®eV < m a < 10 -3 eV and leV < m a < 5eV 
(hadronic axions only). Both these windows are accessible to experiment, and a variety of 
very interesting experiments, all of which involve “heavenly axions,” are being planned or 
are underway. 
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I. Motivation for, and Properties of, the Axion 

Quantum Chromo Dynamics (QCD) is a remarkable theory and is almost universally 
believed to be the theory of the strong interactions. Aside from the calculational difficulties 
associated with actually calculating the spectrum of states in the theory, QCD has but- one 
serious blemish: the strong CP problem. 1 Namely, the fact that non-perturbative effects 
violate CP, T, and P, and unless suppressed would lead to an electric dipole moment for 
the neutron which is in excess of experimental limits by some 10 orders-of-magnitude. 
Proposed in 1977, Peccei-Quinn symmetry with its associated axion is perhaps the most 
elegant solution to this nagging problem. 2 

Non-Abelian gauge theories have a rich vacuum structure owing to the existence of 
non-trivial, vacuum gauge configurations. These degenerate vacuum configurations are 
characterized by distinct homotopy classes that cannot be continuously rotated into one 
another and are classified by the topological winding number n associated with them, 

n= H^f <i3lTr )A’(z )A k (z) 

where g is the gauge coupling, A* is the gauge field, and temporal gauge (A° = 0) has been 
used. The correct vacuum state of the theory is a superposition of all the vacuum states 

I » >, 

|0 >= S„exp(— in0)|n > 

where a priori 0 is an arbitrary parameter in the theory which must be measured. The 
state |0 > is referred to as “the 0- vacuum.” By appropriate means the effects of the 
0-vacuum can be recast into a single, additional non-perturbative term in the QCD La- 
grangian, 

£QCD = £pERT + ®32^2 G * ,lV Gani' 

0 = 0 + Arg det M. 

where G a,iv is the field strength tensor, G afiu is the dual of the field strength tensor 
(GG oc E co \ or • -B co ior ) ) ant ^ ^ is q uar k mass matrix. Note that the effective 0 term in 
the theory involves both the bare 0 term and the phase of the quark mass matrix. Such a 
term in the QCD Lagrangian clearly violates CP, T, and P, and leads to a neutron electric 
dipole moment of order 3 

d n ~ 1O~ 15 0 e — cm 

The present experimental bound to the electric dipole moment of the neutron, 4 d n < 10- 25 
e-cm, constrains 0 to be less than (or of order of) 10 -10 . 
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Before going on to discuss the axion some general comments about the strong-CP 
problem are in order. The unwanted, non-perturbative term in the Lagrangian arises due 
to two separate and independent effects: the 0 structure of the pure QCD vacuum; and 
electroweak effects involving quark masses. In the limit that one or more of the quarks are 
massless the GG term has no physically measurable effects: The 0 term can be rotated 
away by a chiral rotation, and there is no strong CP problem. In the absence of a massless 
quark species (for which the evidence in our world is strong), the effective GG term is 
made of two unrelated contributions which a priori have no reason to cancel. 

One might be tempted to ignore this mysterious topological contribution to the QCD 
Lagrangian, on grounds that one has no need for it, or the hope that its absence will be 
understood at some future date. This is not a particularly good thing to do either; the 
0 structure of the QCD vacuum has at least one beneficial feature: the resolution of the 
U(1)a puzzle. In the absence of such a term one would expect 4 Goldstone bosons when 
the U(2)l ® U(2)r global symmetry (of a massless, 2 flavor QCD world) is spontaneously 
broken by QCD effects. These Goldstone bosons are the 7 r and rj mesons, the rj meson 
being the Goldstone of the spontaneously broken TJ(1)a global symmetry. When non-zero 
u and d quark masses are taken into account one can show that the mass of the 7) must 
satisfy: 5 < y/3m.,r, which, needless to say, is contradicted by reality. The existence of 
the 0 vacuum structure of QCD corrects this erroneous prediction and solves the U(1)a 
problem. 

So it seems likely that the 0 structure of the QCD vacuum is to be taken seriously. 
Moreover, it seems unlikely that any of the quark flavors is massless. How then, is one to 
solve the strong CP problem? The most elegant solution is the one proposed by Peccei and 
Quinn 2 in 1977 (and my personal favorite, as I was a very impressionable young graduate 
student at Stanford when they made their exciting proposal!). Their idea is to make 0 
a dynamical variable, which owing to its classical potential relaxes to zero. This end is 
accomplished by introducing an additional global, chiral symmetry, now known as PQ 
(or Peccei-Quinn) symmetry, which is spontaneously broken at a scale /pq. Weinberg 
and Wilczek realized that because 1^(1 )pq is spontaneously broken there should be a 
Goldstone boson, “the axion” 2 (or as Weinberg referred to it for a while, “the higglet”). 
Because U(1)pq suffers from a chiral anomaly, the axion is not massless, but acquires a 
small mass of order Aq CD //pq, about which we shall be more precise shortly. Moreover, 
due to this anomaly a term in the QCD Lagrangian of the form 

£qcd = - + cons tj— 327[ . 2 G* G ^ 

arises, where a is the axion field and cons' t is a model-dependent constant. Note the 
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similarity of this term to the previously discussed 0 term. These two terms amount to a 
potential for the axion field, which is minimized by having 


< a >= 


Q/pq 

cons't 


for which the coefficient of the offending GG term vanishes! Expanding the axion field 
about its vacuum expectation value < a > the axion part of the QCD Lagrangian is 

2 

£*xion = -d^ad^a/ 2 + cons' t-^-^G^G^ 

(where I have not yet included the axion’s other interactions). Note that the 0 parameter 
has been effectively replaced by the dynamical axion field, whose mass arises due to the 
non-perturbative GG term. 

In an axion model then, the price for resolving the strong CP problem is the existence 
of an additional, spontaneously-broken global symmetry (which often arises in supersym- 
metric and superstring-inspired models in any case) and its associated pseudo-Goldstone 
boson. A priori the mass of the axion (or equivalently the PQ symmetry breaking scale) 
is arbitrary: all values solve the strong CP problem equally well. Taking fpq to be some- 
where between 100 GeV and 10 19 GeV, the associated axion mass then lies between ~ 1 
MeV and 10 -12 eV — a span of some 18 or so orders-of-magnitude to search. 

So much for the high brow theory and philosophy! In order to search for the axion 
one must know about its properties and how it couples to ordinary matter. As eluded to 
above, an axion model has one basic, free parameter: the axion mass, or equivalently the 
PQ symmetry breaking scale. They are related by 

Jl f n m n 0.62 eVIO 7 GeV 

in — : 

1 + z/pq/JV fpq/N 

where z = m u /md ~ 0.56, m x = 135 MeV and / K = 93 MeV are the pion mass and decay 
constant, N is the color anomaly of the PQ symmetry, and our normalization conventions 
for fpQ follow those of refs. 6. The axion field a is related to 0 by: a = (fpq/N)Q. 
The effective Lagrangian for the interactions of axions with ordinary matter (nucleons, 
electrons, and photons) is 

Ant = t|^-^o(!V'7 , 1 75JV) + i^Ld.a(e 7 ^e) + g^aE -B 

with the associated Feynman diagrams shown in Fig. 1. [In all instances where there is only 
one Goldstone boson in the problem, the pseudo-vector coupling may be written instead 
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as a pseudo-scalar coupling; e.g., ig*fixa(N'ysN), by means of a suitable phase rotation of 
the fermion fields. For further discussion of this point, see refs. 7.] 

The axion couplings g * y are given by 


Saee = [XJN + (2a 2 /^)(E\n(f PQ /m e )/N - 1.95ln(A Q c D /m e ))] *— 

(jpq/N) 


9*-n — 


a/2ir 


(. E/N - 1.95) 


g > nn = [(--FUo - F A3 )(X d /2N - 0.18) + (~F A0 + F A3 )(X U /2N - 0.32)]-^- 

jpq/ n 


<7. pp = [(-.FUo - F A3 )(X u /2N - 0.32) + (~F A0 + F A3 )(X d /2N - 0.18)]-^- 

fPQ/N 

where E is the electromagnetic anomaly of the PQ symmetry, 1.95 = 2(4 + z)/ 3(1 + z), 
0.32 = 1/2(1 + z), 0.18 = z/ 2(1 + z), a ~ 1/137 is the fine structure constant, F A o — 
—0.75 is the axial- vector, isoscalar part of the pion-nucleon coupling, and F A3 ~ —1.25 
is the SLxial-vector, isovector part of the pion-nucleon coupling. [Note that the axion- 
nucleon couplings have been computed in the context of the “naive quark model;” more 
sophisticated treatments lead to slightly different axion-nucleon couplings (see, e.g., Mayle, 
et al. 35 ).] 

The quantities Xi (t = u, d, e ) are the PQ charges of the u and d quarks and the 
electron. Depending upon the PQ charge of the electron, the axion can be classified as 
one of two generic types: hadronic, 8 the case where X t = 0 (no tree level coupling to the 
electron); or DFSZ,® the case where all the Xi axe of order unity. [In the original DFS 
model, 9 N = 6, X e /N = cos 2 /3/3, X u /N = 1 — cos 2/?, and X d = 1 + cos 2(3. Here 
parameterizes the ratio of the ‘up’ and ‘down’ PQ vacuum expectation values.] 

First note that all the axion couplings (flf*u) are proportional to 1 /(/pq/N), or equiv- 
alently m, : the smaller the axion mass, or the larger the PQ SSB scale, the weaker the 
axion couples. The coupling of the axion to the photon arises through the electromag- 
netic anomaly of the PQ symmetry, and allows the axion to decay to two photons, with a 
lifetime, 


T a = 6.8 x 10 24 sec; 


(m./eV) 


-5 


[(E/N — 1.95) /0.72] 2 

Note too that the coupling of the axion to 2 photons depends upon the ratio of the 
electromagnetic to the color anomaly; when the axion is incorporated into the simplest of 
GUTs, E/N = 8/3 and (E/N — 1.95) ~ 0.72. However, it is possible that E/N could have 
a different value, even 2, in which case its 2 photon coupling (E/N — 1.95 ~ 0.05) would be 
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strongly suppressed. 10 We should keep this fact in mind, as it will be of some importance 
when discussing the astrophysics! effects of a hadronic axion. 

Next note that the coupling of the axion to the electron has a tree level contribution 
which is proportional to X e and vanishes for the hadronic axion, and a loop correction 
which is proportions! to a 2 and arises due to the anomalous 2 photon coupling of the 
axion (see Pig. 1). The axion-electron coupling is of great importance in determining the 
astrophysical effects of the axion, and explains why it is usually necessary to discuss the 
astrophysics! constraints to DFSZ and hadronic axions separately. 

Finally, note that the axion-nucleon coupling arises from two roughly equal contribu- 
tions: the tree level coupling of the axion to up and down quarks, and a contribution 
which arises due to axion-pion mixing (both the axion and pion are Goldstone bosons with 
the same quantum numbers, and the physically-propagating states mix). This means that 
even a hadronic axion which does not couple to light quarks at tree level (as was the case 
with the original hadronic 8 axion which only coupled to one very heavy, exotic quark) still 
has a coupling to nucleons which is comparable to that of a DFSZ axion. Because of this 
fact, the bound to the axion mass based upon SN 1987A, which involves the axion-nucleon 
coupling, is essentially the same for the hadronic and DFSZ axion. 

While I have taken the view here that a priori the axion mass (or symmetry breaking 
scale) is an arbitrary parameter, to be determined by experiment, that viewpoint belies 
the history of the axion. The original axion proposed by Peccei and Quinn was based upon 
a PQ symmetry breaking scale equal to that of the weak scale (fpQ ~ 250 GeV), leading 
to an axion mass of about 200 keV. As I will now discuss, such an axion was quickly 
r\iled out by unsuccessful experimental searches. Shortly thereafter the ‘invisible axion’ 
was invented, 8 ’ 8 an axion with symmetry breaking scale 250 GeV and mass <C 200 keV, 
whose interactions are necessarily extremely weak (recall g <x m a <x f pq)’ And of course, 
it goes without saying that once the weak scale had been ruled out, all educated bets as 
to the PQ SSB scale were off! 

[Very recently it has been argued that wormhole effects, might cause the wave function 
of the Universe to be very highly peaked at 0 = 7r (a CP-conserving value); if correct, this 
would solve the strong-CP problem and obviate the need for the axion. 11 Even if this is 
true, Nature may still provide us with an axion, as PQ symmetry seems to be very generic 
to supersymmetric and superstring inspired models.] 
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II. Laboratory Searches 


As mentioned above, the original Peccei-Quixm- Weinberg- Wilczek axion 2 was charac- 
terized by a SSB scale of order the weak scale and a mass of order 200 keV. As such its 
interactions were roughly semi-weak, making it accessible to laboratory searches. So ac- 
cessible in fact, that the original axion was very quickly ruled out. Without doing justice 
to the history and to the variety of important experiments, I will very briefly mention the 
most sensitive laboratory searches. 12 First is the Kaon decay process 

K+ _> 7T+ + a 

where the axion goes unseen. The present experimental upper limit to the branching ratio 
for K* — » 7r + + nothing is 3.8 x 10 -8 . 13 In an axion model this process arises either through 
axion/pion mixing (the decay K + — ► ?r + + 7r° is observed), or an off diagonal coupling of 
the axion to s and d quarks. 

Next are the decays of quarkonium {QQ) states 

J/ij> a + 7 T — >0 + 7 

The upper limit for the branching ratio for these two processes are 1.4 x 10~ 5 and 3 x 10“ 4 
respectively. 14 

Based upon the three processes just mentioned, one can safely conclude that 

fpQ > 10 3 GeV or m 0 < 6 keV 

There have been other axion searches involving disallowed J p = 0“, 1 + , 2 ~, ... nuclear 
transitions, reactor experiments, and beam dump experiments. In general, the limits that 
follow from these are less strigent or more difficult to interpret. 

While laboratory-based experiments exclude axion masses in the range of about 10 
keV to 1 MeV and most certainly rule out the original axion, they leave open an enor- 
mous window: lO -12 eV — lOkeV, one which has only been explored by astrophysics and 
cosmology. 


III. Axions and Stars 

The life of a 6tar is rather uncomplicated, the simple struggle to lose the enormous 
nuclear free energy associated with the primordial composition of the Universe (for every 
10 atoms, roughly 9 H atoms, 1 He atom, and a trace of D, 3 He, and 7 Li). Given the 
intrinsically short time scale associated with nuclear interactions, it should be a great 
surprise that stars live as long as they do: a star like our sun bums hydrogen for 10 10 yrs. 
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The reason for the hang up is simple; the rate at which a star can liberate its nuclear 
free energy is controlled not by nuclear reaction rates, but rather by the rate at which the 
nuclear energy liberated can be transported through the star and radiated into the vacuum 
of space. Under the conditions that exist in a typical star, say our sun, the mean free path 
of a photon is only about a cm!, and the time required for a photon liberated at the center 
of the sun to make its way (figuratively of course) to the surface is of order 10 7 yrs. The 
enormous opacity of ordinary matter to photons of course traces to the strength of the 
electromagnetic interactions; recall that the Thomson cross section is a ~ 0.67 x 10~ 24 
cm 2 . The very long time required for a star like our sun to burn its hydrogen fuel then 
owes to the large interaction cross section of the photon with ordinary matter. 15 

The existence of a light (i.e., compared to typical stellar temperatures, T ~ keV-MeV), 
weakly-interacting particle has the potential to greatly accelerate the evolutionary process 
of stars of all types by more efficiently transporting energy away, and thereby to shorten 
their lifetimes. To effectively carry off the free energy liberated in the nuclear reactions in a 
star, the hypothetical “super coolant” must interact weakly enough so that it streams right 
out without interacting, but strongly enough so that it is produced in sufficient numbers to 
carry away large amounts of energy. As one might guess, the optimal interaction strength 
is such that the super coolant particle has about 1 interaction as it streams out. Nature 
has provided us with at least 3 candidates, the 3 neutrino species, and comptemporary 
theorists have postulated another, the axion. 

Before turning to the axion, let us orient ourselves by discussing neutrino cooling in 
stars. Because of the nature of the weak interaction neutrino cross sections are highly tem- 
perature sensitive, proportional to GpT 2 . [It is interesting to note that neutrino emission, 
unli ke axion emission, is necessarily a second order weak process, i.e., C„ oc< v > -4 , where 
< V >~ 250 GeV is the SSB scale of the weak interactions.] In ordinary main sequence 
stars the neutrino luminosity C v is proportional to T®, whereas we shall see shortly that 
the photon luminosity jC t oc TV 2 (T c is the central temperature of the star). Only in stars 
hotter than about 10® K does neutrino cooling begin to compete with photon cooling; for 
these stars (O, Si burning stars and beyond) neutrino emission is the dominant cooling 
mechanism, and as a result the time scale for these burning phases is greatly reduced as 
neutrinos can just stream out (O burning time scale is of order 10 5 yrs; Si burning time 
scale is of order sec’s). In fact, long ago it was argued that the existence of carbon burning 
stars places a limit to Gjr; 16 had Gp been a factor of 3 or so larger, the evolution time scale 
for C stars would have been greatly reduced due to neutrino emission, so much so that C 
burning stars would evolve through C burning so quickly that none would be observed. 
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As a preview of things to come, we should mention SN 1987A; the primary cooling 
mechanism for the hot, nascent neutron star was, as we all now appreciate, neutrino 
emission. Moreover, based upon the neutrino burst time scale, it has been argued that 
the number of light neutrino flavors (mass < MeV) must be less than about 9, otherwise 
neutrino cooling of the neutron star would have proceeded more rapidly (by a factor of 
roughly 3) than observed. 28 

The effect of axion emission on stars is clear: the acceleration of their evolution and 
shortening of their lifetimes. In main sequence stars and red giants the primary axion emis- 
sion processes are: 17 the Compton-like process 7 + e - — ► a + e - ; and axion bremsstrahlung 
e~ + Z—*a + e~+Z t both of which are proportional to g* ee oc m 2 . Of lesser importance 
unless g aee vanishes at tree level, as it does for a hadronic axion, is the Primakoff process 
7 + Z( or e“) — ► a + Z( or e”). In very low-mass stars (M < 0.2 M@) emission through 
the axio-electric effect (the analogue of the photo-electric effect) is also very important. 18 

To begin, consider the sun. At the center of the sun the temperature is about 1.6 x 10 7 
K, and nuclear reactions (p + p — > D + e + + v t \ D + D —* 4 He + 7) liberate free energy at 
the rate of a few ergs g -1 sec -1 . On the other hand, axion emission carries away energy 
at the rate of 

£„ ~ lerg g -1 sec -1 (T e /10 7 K) 8 [10 7 GeV/(/pQ/JV)] 2 

Roughly speaking then, if (/ pq/N ) were less than about 10 7 GeV, axions would carry 
energy away from the center of the sun faster than nuclear reactions could generate it. 
[Note the above rate is that for a DFSZ axion.] The thermal time constant of a star 
like the sun is only about 10 7 yrs; this is the time required for the star to radiate away 
its thermal energy reserves, and is known as the Kelvin- Helmholtz time. Thus an axion 
luminosity greater than the rate at which nuclear energy is released can only be tolerated 
on a short time scale (< 10 7 yrs, or so); if (/pq/N) were less than 10 7 GeV then the 
sun would have to “adjust itself” to re-establish energy balance. As we shall see, in a 
hypothetical star in which one “turns on” axion emission, the star contracts to raise its 
temperature and nuclear energy liberation rate to balance axion losses. In the process it 
would also raise its photon luminosity, and as a result of both axion emission and enhanced 
photon emission its lifetime would be shortened. Thus, the all important observable for 
constraining axion emission from the sun is its age at a given 4 He mass fraction. For 
( fpQ /N) < 10 7 GeV, a sun with our sun’s 4 He abundance would be younger than our sun 
is known to be. 19 

Let us consider the sun in slightly more detail. As we have discussed the photon 
luminosity of the sun is determined by the opacity of solar material. Just by analyzing hy- 
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drostatic equilibrium and energy transport in stars like the sun (i.e., stars less massive than 
about 2 Mo), Chandrasekhar has derived a remarkable formula (the so-called luminosity 
formula) which relates the photon luminosity of a star to its central temperature 15 

£ 7 oc ( GfjL e ) 7 M*T t 'l> 2 

where /z e is the mean molecular weight per electron (1 for a pure H star; 2 for a pure 4 He 
star), M is the mass of the star, and T c is the central temperature. 

Energy balance requires that the energy liberated by nuclear reactions 

Qnuc = / inucdM 
Jt tar 

be equal to the photon (plus axion) luminosity. In general the nuclear energy liberation rate 
(here per g of material per sec) is very temperature dependent and can be parameterized 
as 

Snuc « pT n oc T n+3 

where for the sun n ~ 3 (in a star like the sun the entropy per baryon T 3 /p is constant, 
so that p oc T 3 ). To begin, consider a star in the absence of axion emission. Energy 
equilibrium requires that £ 7 = Q nuc = Qo* Now suppose that the star radiates axions, 
with an axion luminosity Q a = cQ 0 . Energy balance now implies that 

C-f + Q a = Qnuc 

6£-y Q a = &Qnuc 

By fiat Q a = cQo, and using Chandrasekhar’s luminosity formula we see that $£ 7 = 
0.5(6T e /T c )Qo • Finally, since e nuc oc T n+3 , we have ^Q„ U c = (n + 3)(£T c /T c )Qo- Using the 
perturbed energy balance equation we find that 

6T C /T C = e/6.5 5£ 7 /Q 0 = c/13 

6Ctot/Q 0 = 14e/13 SR/R 0 = -6T C /T C = -e/6.5 

where we have used the fact that the radius of the sun R oc T~ x to find SR/R 0 . 

As advertised, we see that a star perturbed by axion emission contracts to raise its 
temperature and restore energy equilibrium, and in the process increases its photon lu- 
minosity also. Suppose that e = 1/2; then the central temperature increases by about 
8%, and the total luminosity by about 54%, thereby decreasing the evolution time scale 
by more than a factor of 2, strongly suggesting that e > 1/2, or so, is inconsistent with 
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our knowledge of the sun. [To make such an argument rigorous one must also consider 
the very strong compositional dependence of stellar models, reflected in £ 7 a /z 7 ; this has 
been done in ref. 19.] 

An even more sensitive barometer for stellar axion emission is the rate of 8 B neutrino 
emission 20 (the high energy neutrinos which have been detected by Davis’ 37 C1 experiment 
and by the KII detector). The rate of emission of these high-energy neutrinos is propor- 
tional to T £ ( p ~ 13, depending upon which quantities in the stellar model are held fixed). 
Using our previous formula for ST C we see that the 8 B neutrino flux would increase by 
almost a factor of 3 (for e = 1/2), exacerbating an already large discrepancy. We see that 
the simple-minded limit provided by i a < e nuc is more than justified! 

Finally, it is interesting to note that Germanium double beta decay experiment of 
Avignone, et al. 21 provides a similar limit to / pq/N based upon the non-observation of 
solar axions in their detector. [Solar axions would be detected by their interactions with 
electrons in the Ge detector: e~ + a —* e~ + 7.] 

To summarize the axion mass limits that follow from the sun, they are: m 0 < 1 eV 
(DFSZ); and m 0 < 20/[(.E/jV — 1.95)/0.72] (hadronic) (see Fig. 2). Moreover, these limits 
do not apply to an axion of mass greater than 10 keV or so, since the production of such 
axions would be severely suppressed owing to the fact that the temperature at the center 
of the sun is only a few keV. While we have explicitly displayed the dependence of the 
hadronic axion’s model dependent coupling to two photons, we have not been so careful 
with the model dependence of the DFSZ axion’s coupling to electrons. DFSZ axion mass 
limits are necessarily proportional to cos -2 0, and in the limit that (3 — » n/2 ( X e — ► 0), 
they revert to those of the hadronic axion. 

The discussion above should provide the reader with the flavor of stellar limits to the 
axion mass. They all rely on the fact that axion emission modifies stellar evolution in such 
a way as to significantly affect an observable, usually the lifetime of the star. Now let us 
turn to the most stringent stellar evolution limits that exist at present. These limits are 
provided by the evolution of red giant stars, 6tars whose centred temperatures reach 10 8 
K and whose central densities are ~ 10 2 — 10 4 gcm -3 . Because hadronic and DFSZ axions 
couple very differently to electrons, their mass limits are very different. 

The constraint to the mass of the hadronic axion is based upon the helium-burning 
lifetimes of red giant stars. 24 As we discussed above, when axion emission is taken into 
account the central temperature of the star is necessarily increased to satisfy the extra 
energy being carried away by axions, and this accelerates the evolution and shortens the 
lifetime. The helium- burning phase of a red giant lasts of order 10 8 yrs or so — too long for 
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most astronomers to observe. However, when one observes a cluster of stars (say M67, for 
example), the number of helium-burning red giants one sees is determined by the length 
of time red giants typically spend burning helium — the shorter the time, the fewer that 
will be seen. Raffelt and Dearborn 22 argue that a hadronic axion of mass greater than 
about 2eV /[(E/N — 1 .95) /0.T2] would reduce the helium-burning time scale by more than 
an order of magnitude, in severe contradiction with observations of the number of helium- 
burning red giants seen in the cluster M67. Once again we see the factor which arises from 
the model dependence of the axion-photon-photon coupling. 

The red giant limit for the DFSZ is based upon a slightly more subtle dynamical 
argument. 23 Before helium ignition occurs in the core of a red giant, the 4 He core is 
supported by electron degeneracy pressure. This is a very dangerous condition because 
any increase in temperature is not accompanied by a similar increase in pressure, and so 
once any 4 He is ignited nuclear burning is a runaway process, until thermal pressure support 
becomes dominant. The brief period of thermal runaway is referred to as the helium flash 
(not to be confused with a hot flash). [In an ordinary star, the simple physical fact that the 
pressure is proportional to the temperature stabilizes nuclear burning, as any increase in 
temperature is accompanied by an increase in pressure which causes the star to expand and 
thereby cool — Nature’s stellar thermostat!] Before the helium flash, hydrogen continues to 
burn just outside the helium core. As the helium core grows in mass, its radius decreases 
(for degenerate matter R oc M -1 / 3 ), and the accompanying release of gravitational binding 
energy heats the core. Eventually, the helium core becomes hot enough for the triple-a 
process to burn helium to carbon. The effect of axion cooling in the helium core decreases 
the temperature rise in the helium core associated with the contraction, and according to 
Dearborn, et al. 23 , can prevent helium ignition from ever taking place for a DFSZ axion 
of mass greater than about 10 -2 eV. 

Let us consider their argument in slightly more detail. As the helium core increases in 
mass the size of the core contracts ( MR 3 = cons't), thereby releasing gravitational energy: 



a M l **M 


In the absence of axion cooling, the dominant cooling mechanism for the approximately 
isothermal core is neutrino emission (because of the long mean free paths of electrons in 
degenerate matter, degenerate matter is almost always isothermal). For our purposes, let 
us assume that axion emission dominates, which is the case for an axion of mass sufficient 
to “screw up” the helium flash. The Compton-like process is the dominant axion emission 
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process in the helium core, and the energy radiated in axions is 

Q a ~ mlMT* 

Energy equilibrium, i.e., E grav = Qa, determines the temperature of the core: 

Tcore ~ 

Note, the larger the axion mass the lower the temperature of the core; this differs from 
the usual case where the existence of axions actually causes the star to raise its temperature 
to compensate for axion emission. Of course, in the present situation that is not possible 
since nuclear reactions are not yet occurring in the 4 He core. One can easily appreciate how 
a sufficiently massive axion can prevent the core from reaching the temperature required 
to ignite helium. However, the unstable nature of degenerate matter gives one pause; to 
ignite the core one only has to “light” the smallest region, which then triggers thermal 
runaway and ignites the entire core. Since the SN 1987A bound will be more stringent 
than this one, we do not have to lose sleep over this point. 

To summarize the red giant constraints, the helium-burning lifetime argument pre- 
cludes a hadronic axion of mass greater than about 2 eV/[(E/N — 1.95)/0.72], while the 
ignition of helium burning precludes a DFSZ axion more massive than about 10~ 2 eV (see 
Fig. 2). Because the temperatures in the cores of red giant stars are of order 10 keV or so, 
emission of axions more massive than about 100 keV is severely suppressed, and so these 
arguments do not apply to axions more massive than about 200 keV (which of course are 
precluded by laboratory searches). 

Before going on to discuss the most stringent astrophysics! bound, that based upon the 
cooling of yet another kind of star, the newly-born, hot neutron star associated with SN 
1987A, we should mention two other astrophysical bounds. First, the cooling of relatively 
young (few 100 yrs old) neutron stars, including the Crab pulsar and RCW 103. It has been 
argued that an axion (of either type) more massive than about 10 -2 eV or so would cool 
several of the well known young neutron stars so rapidly as to be inconsistent with Einstein 
measurements of their surface temperatures. 24 Similar arguments based upon the cooling 
of white dwarf stars seem to preclude a DFSZ axion of mass greater them about 3 x 10“ 2 
eV. 25 However, for both of these bounds there are theoretical, as well as observational, 
uncertainties which cast some doubt upon them. In the case of the neutron star bound, 
it is not clear whether thermal emission has actually been detected from all, or even any, 
of these young neutron stars. On the theoretical side, it has been noted that superfluidity 
in cool neutron stars would greatly reduce axion emission and significantly degrade the 


13 


constraint. There are similar residual theoretical and observational uncertainties associated 
with the white dwarf bound. Since other constraints with fewer uncertainties exist that 
are as strong, or stronger, we will not dwell on the reliability of these limits further. 

The various astrophysical constraints to the axion mass based upon stellar evolution 
are summarized in Table 1. 


IV. Axions and SN 1087A 

SN 1987A not only confirmed astrophysicists’ more cherished beliefs about type II (core 
collapse) supernovae, but also provided a unique laboratory for the study of the properties 
of ordinary neutrinos, right-handed neutrinos, axions, and other exotic hypothetical par- 
ticles. Here we will be interested in the bound to the axion mass provided by SN 1987A. 
As we shall see, it is the most stringent, and I believe, the most reliable astrophysical 
constraint to the axion mass. 

Shortly after the gravitational collapse and hydrodynamic bounce of the 1.4 M© Fe 
core of the blue super giant Sanduleak -69 202 (thought to be a ~ 15 M© stax), the 
central temperature of the nascent neutron star was 20-70 MeV and the central density 
was ~ 8 x 10 14 g cm -3 . During the catastrophic collapse of the Fe core about 3 x 10 53 
ergs of binding energy were liberated, and according to the standard picture, this energy 
is radiated in thermal neutrinos of all 3 types. 38 The neutrino mean free path within the 
core is much smaller than the size of the core (~ 10 km) and so even neutrinos are trapped 
in the core. Thus, neutrinos are radiated from a neutrino sphere ( R ~ 15 km, p ~ 10 12 g 
cm -3 , T ~ 4 MeV). 

Neutrino emission is characterized by two phases: the first is powered by residual 
accretion and hydrodynamic contraction of the outer core, and lasts 1-2 sec; the second 
phase is powered by the diffusion of heat trapped in the inner core region, and lasts ~ 5 — 10 
sec, the time scale for neutrino diffusion from the core to the neutrino sphere. The energies 
associated with the two phases are comparable, and as a result one expects a neutrino 
burst of order 5-10 sec. The observations of KII 27 and IMB 28 are both qualitatively and 
quantitatively consistent with the standard picture. 28 As it does with other types of stars, 
if the axion exists, it can play an important role in the cooling of this nascent neutron stax. 
In the case of SN 1987A the observable effect of axion cooling would be the shortening of 
the neutrino burst, and fortunately we have 19 beautiful neutrino events spread out over 
about 10 sec to study the potential effects of axion cooling! 

Under the conditions that existed in the post collapse core the dominant axion emission 
process is nucleon-nucleon, axion bremsstrahlung. 29 In the one pion exchange approxima- 
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tion (OPE) there are 4 direct and 4 exchange diagrams (see Fig. 1). The relevant axion 
coupling here is that to nucleons, which we should we recall is relatively insensitive to the 
type of axion and is of order tun l{fa/N) ~ 10 -7 (m a /eV). The full matrix element squared 
(64 terms) has been evaluated exactly in the OPE approximation. 30 From |A4| 2 the axion 
emission rate (here, per volume per time) is given as 

i a = J dn 1 dn 2 dn3 < fn 4 dn 0 (27r) 4 5|M| 2 5 4 (p 1 + P2 - j> 3 - pa - p«)£«/i/ 2 (i - /*)(i - / 4 ) 

where dllj = d 3 pi/(2ir) 3 2Ei, the labels i = 1—4 refer to the incoming (1,2) and outgoing 
(3,4) nucleons, i = a denotes the axion, 5 is the symmetry factor for identical particles 
in the initial and final states, |A4| 2 is summed over initial and final nucleon spins, and 
the nucleon phase space distribution functions are fi = [exp (Ei/T — m/T) + l] -1 . The 
emission rate is relatively easy to evaluate in the fully degenerate or non degenerate regimes; 
however, the nucleons in the core are semi-degenerate, cfermi ~ T. In addition, since the 
post collapse core has roughly equal numbers of neutrons and protons, 3 bremsstrahlung 
processes are important: nn — ► nn + a, pp —* pp + a, and np —* np + a. The axion 
emission rate, for all 3 processes and arbitrary nucleon degeneracy, has been evaluated 
numerically. 30 [As it turns out, the non- degenerate axion emission rate provides a good 
approximation to the actual rate for the semi-degenerate conditions that exist in the hot 
neutron 6tar. 30 ] 

Axions less massive than about 0.02 eV, once radiated, freely stream out of the nascent 
neutron star, and thereby accelerate the cooling. Qualitatively then, one would expect 
axion emission to shorten the duration of the neutrino pulse — this is in fact what occurs 
(see Fig. 3). [Of course axion emission, which proceeds predominantly from the high 
temperature, high density inner core does not directly affect neutrino emission, which 
proceeds from the neutrino sphere (in the outer core).] We have incorporated axion cooling 
into realistic numerical models of the initial cooling of the nascent neutron star; 31 the 
biggest theoretical uncertainty in these models is the equation of state (EOS) above nuclear 
density, densities which are achieved in the core during and after collapse. We have allowed 
for a wide range of EOS’s, from a very stiff EOS to a very soft EOS. For our various axion- 
cooled, numerical models we have computed the resulting neutrino flux and the predicted 
response of the KII and IMB detectors: expected number of events; and burst duration, 
At(90%), the time required for the number of events to achieve 90% of its final value. The 
quantity Al(90%) is the most sensitive indicator of axion emission. Axion emission tends 
to rapidly cool the inner core, depleting the energy which powers the second part of the 
burst. This effect is clearly seen in Fig. 3 where A/(90%) is plotted as a function of m 0 . 
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Axion emission has virtually no effect on At(90%) until a mass of ~ 3 x 10~ 4 eV, and by 
an axion mass of 10 -2 eV the duration of the neutrino burst has dropped to less than a 
sec (~ time scale for the first phase of the burst). For comparison, for an axion mass of 
10 -2 eV, the expected number of neutrino events has only dropped from ~ 10 to ~ 8 for 
KII and from ~ 6 to ~ 4 for IMB. Likewise, for an axion mass of 10~ 2 eV, neutrinos still 
carry away more than 50% of the binding energy. The large effect on the burst duration 
traces to the fact that axion emission from the core efficiently radiates away the heat which 
powers the latter phase of the burst. 

[One might wonder if a finite mass for the electron neutrino could lengthen an axion- 
shortened neutrino burst. A 20 eV or so mass might work just fine for the KII events; 
however, because the energies of the IMB events are much larger on average, a mass of 
30-50 eV would be required to lengthen the IMB burst, a value precluded by the KII data 
and laboratory experiments.] 

For axion masses greater than ~0.02 eV axions interact sufficiently strongly so that 
they do not simply stream out: rather, they become trapped in the core and are radiated 
from an axion sphere, with temperature T a . In the trapping regime, the axion luminosity 
is oc T*. With increasing axion mass, the axion sphere moves outward and therefore 
has a lower temperature. Thus, for m a > 0.02 eV the axion luminosity decreases with 
increasing axion mass; whereas in the freestreaming regime the axion luminosity increases 
with increasing axion mass, as m 2 (see Fig. 4). For sufficiently large axion mass the effect 
of axion cooling becomes acceptable. The complexity of axion transport has thus far 
prevented us from incorporating axion cooling into our numerical models (although such 
work is in progress). However, simple analytical models indicate that for an axion mass 
of ~ 2 eV or greater the axions are so strongly trapped that their presence is equivalent 
to less than a couple of additional neutrino species and is therefore consistent with the 
observations of KII and IMB. 32 

Two uncertainties cast a shadow of doubt on this limit: the equation of state at super- 
nuclear density and the calculation of the axion emission rate. While the former is indeed 
an important uncertainty, we have explored a variety of EOS’s and our limit does not vary 
significantly. The latter is of greater worry. The axion emission rate has been calculated 
in the OPE approximation at supernudear densities, neglecting any finite density effects. 
Broadly speaking then, there are two concerns: the validity of the OPE approximation it- 
self; and possible effects due to finite density. The validity of the OPE approximation has 
been addressed by calculating the cross section for the analogue process, pp — ► pp + 7r°, us- 
ing OPE, and comparing the result to the existing experimental data. 33 ’ 34 The agreement 
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between the OPE cross section and the experimental data is quite impressive (better than 
a factor of 2) at the energies of interest. Since both the axion and pion are derivatively- 
coupled Goldstone bosons, one would expect this analogue process to provide a good test 
of OPE for axion production, a test which OPE passes with flying colors. 

The finite density effects are more difficult to access. However, using the non-linear 
sigma model as guide for the behaviour of the interactions of pions and nucleons at high 
density, the authors of ref. 34 have concluded that any such effects are likely to be small 
(less than a factor of 3), and could possibly enhance axion emission, thereby strengthening 
the SN 1987A axion mass bound. The one finite density effect which might significantly 
modify the SN 1987A bound is the existence of an exotic state of matter at the core of the 
neutron star; e.g., a pion condensate, quark matter, or strange matter. Although such a 
possiblity seems to be a long shot, we are currently studying how the axion limit would be 
affected. 

In sum, the SN 1987 A axion mass constraint is the most stringent astrophysical con- 
straint, and applies equally to both types of axions. The inadequacies of OPE and the 
effects of finite density together might account for a factor of 3 or so uncertainty in the 
axion emission rate. Since that rate itself scales as the axion mass squared, that uncer- 
tainty probably amounts to less than a factor of y/3 uncertainty in the axion mass limit 
which follows. Moreover, because the physics of the cooling of the nascent neutron star 
is so simple and the observable (the neutrino burst) is so clean and direct, the SN 1987A 
constraint is probably the most reliable astrophysical bound. 

[While I have only discussed the work I have been involved in with regard to axions 
and SN 1987A, si mil ar work has been carried out by other authors, 35 and at present all 
axe in agreement as to the mass constraint which pertains in the free streaming limit, i.e., 
m 0 < 10 -3 eV. Thus far the other authors have not addressed the trapped regime.] 

V. Axions and Cosmology 

The topic of relic particles from the early Universe is a most interesting one. When the 
relic being considered is the axion, the topic is even more interesting! Relic axions arise 
due to three, different and distinct processes: thermal production; 38 coherent production 
due to the initial misalignment of the axion field; 37 and the decay of axionic strings. 38 Each 
of these three processes can be the dominant production mechanism, depending upon the 
axion mass and whether or not the Universe ever underwent inflation (see Fig. 5). Let us 
consider the three processes in turn. 

Owing to the state of near thermal equilibrium that existed in the early Universe 
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all kinds of interesting particles were present in great abundance at early times. Roughly 
speaking the criterion for a particle species to be in thermal equilibrium is that the reaction 
rate T for processes which create and destroy that particle species occur rapidly compared 
to the expansion rate of the Universe, H ~ T 2 /m p \, i.e., T > H. For axions the important 
creation and destruction processes are: photoproduction 7 + Q — ► (? + a (Q is a heavy 
quark), Primakoff production 7 + q — * q + a (q is any charged particle), and nucleon- 
nucleon, axion bremsstrahlung N + N—^N + N + a ( N is a nucleon) (see Fig. 1). Since 
each of these processes involves a single axion coupling, the rate for any of them scales as 
m a . 

Based upon a careful analysis of the rates for these processes we find that for an 
axion more massive than about 10 -3 eV there is a period in the history of the early 
Universe, from a temperature greater than a few GeV (or even more for heavier axions) 
down to 100 MeV or so, where axions were in good thermal contact with the universal 
plasma, and should have been present in numbers comparable to photons (more precisely, 
n a EQ = n 7 /2). 38 ’ 39 Then, just as light neutrinos do, thermal axions decouple while they 
are still very relativistic, and their abundance per comoving volume “freezes in”. Today, 
they should have an abundance relative to photons of order n 0 ~ 2n 7 /ll ~ 100 cm -3 . 
[More precisely, n a / s = 0.278/^, (Tj) where Td is the temperature at which axions finally 
decouple and s is the entropy density, today s ~ 7 n 7 ~ 2970cm -3 .] 

The thermal axion contribution to the present mass density of the Universe scales as 
m 0 and is given by 

flthermalh 2 ~ 10 -2 (m a /eV) 

where lOOhkmsec -1 Mpc -1 is the present value of the expansion rate (see Fig. 5). We see 
that thermal neutrinos can only close the Universe for an axion mass of order IOO/ 1 2 eV, 
for which the axion lifetime is shorter than the age of the Universe. [It is intriguing to 
note that for h < 1/2 and E/N ~ 2, thermal hadronic axions would close the Universe, 
would not have decayed by the present epoch, and would escape the previously discussed 
astrophysical bounds.] 

While thermal axions probably cannot provide a significant fraction of the present 
density, it has been pointed out that the lifetimes of multi-eV axions are “well-matched” 
to the present age of the Universe, i.e., sufficiently long so that not all the relic axions have 
decayed, and sufficiently short so that a substantial fraction tire decaying at present. 40 
Since the decay of an axion is a 2-body process, the decay-produced photons are mono- 
energetic (but slightly-broadened due to any velocity that the decaying axions may have, 
A A ~ (v/c)A). While multi-eV thermal axions will not contribute substantially to the 
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present mass density of the Universe, they will find their way into the many gravita- 
tional potential wells that exist, e.g., in galaxies (including our own) and in clusters of 
galaxies. 40 ’ 41 Moreover, their decays will produce potentially detectable line radiation, at 
a wavelength A a = 2 hc/m a = 24800A/(m o /eV). The intensity of this radiation from the 
halo of our own galaxy should be 41 

/halo = 2 x 10 -23 ergcm -2 sec -1 arcsec -2 A -1 (m 0 /eV) 10 [(E/N — 1.95)/0.72] 2 J(8) 

where J{8) is the angular dependence of the signal which owes to the fact that we do not 
reside at the center of our galaxy. The width of such a line is expected to be of order the 
virial velocity in our galaxy, or AA ~ 10 -3 A a . When one turns a telescope to the blank 
night sky one sees many lines — not all from axion decays, rather from airglow! Shown in 
Fig. 6 is a high-resolution spectrum of the night sky (air glow) 42 and the expected axion 
line (for ( E/N — 1.95) = 0.72); from the existing data it is clear that an axion mass of 
greater than 4 eV or so is definitely precluded 41 — so much for thermal axions providing 
closure density. 

Even more favorable is the axion-produced line from relic thermal axions which reside 
in clusters 41 

/cluster = 2 x 10 -20 erg cm -2 sec -1 arcsec -2 A -1 (m a /eV) 7 [{E/N — 1.95)/0.72] 2 

Here the line width is expected to be of order the virial velocity in a cluster, or AA ~ 
10 -2 A o . Moreover, by observing a cluster one has two other advantages: first, one can 
remove many of the airglow lines by subtracting “off cluster” measurements from “on 
cluster” measurements; second, the wavelength of the cluster axion line depends upon 
the red shift of the cluster, A a (cluster) = (1 + Cluster )A a , and by looking at two or more 
clusters with different red shifts one can further discriminate against other night sky lines. 
Currently an observational effort is being mounted to search for a cluster axion line, and 
it is hoped that this effort will be sensitive to an axion as light as about 2 eV. [The 
rapidly decreasing strength of the axion line with decreasing axion mass, together with the 
increasing glow of the night sky at longer wavelengths, precludes searching for axion line 
radiation for masses smaller than about 2 eV.] 

Axions are also produced by a very interesting and highly non-thermal process involving 
the relaxation of the 0 angle. Today the axion mass anchors 0 at the CP-conserving 
value 0 = 0. The axion mass, which arises due to instanton effects, is very temperature 
dependent, 43 

m a (T) ~ 0.1m a (Aqcd/Z’) 3 ' 7 
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At very high temperatures the axion mass is essentiaUy zero. Specifically, at T ~ fpQ, 
when PQ symmetry is spontaneously broken, the axion mass is for all purposes negligible 
and the axion is a Goldstone boson. That means that no special value of 0 is specified 
by dynamics and all values of 0 are equally palatiable! Therefore, the initial value of 0 
must be choosen by some stochastic process, and in general the initial value of 0, call it 
d \ , is likely to be of order unity. Thus at early times the axion field is misaligned with the 
minimum of its potential (0 = 0). 

When the axion mass does “turn on” and become comparable to the expansion rate 
of the Universe, the axion field will start to roll toward 0 = 0, and of course will over- 
shoot 0 = 0. Thereafter, it will oscillate like as heavenly harmonic oscillator. These 
cosmic oscillations of the axion field correspond to a zero momentum condensate of ax- 
ions (with phase space density well in excess of 10 30 )! It is simple to estimate their 
present mass density. The initial energy density trapped in the misaligned axion field 
is, p a = m 0 (Ti) 2 aJ = to o (Ti) 2 0 2 (/pq/W) 2 . The initial axion number density is just 
Pa/m B (Ti), or 

».(T ^-maiT^OKfpQ/N) 2 


where T\ is the temperature when the axion field begins to oscillate, m a (T \ ) ~ 3H(T \ ) ~ 
T\/m p i. For an axion of mass 10” s eV, T\ is of order a GeV, and T\ scales as 
Assuming that there has been no entropy production since the axion field began to oscillate, 
the axion number density to entropy density ratio is conserved (even in the presence of the 
time- varying axion mass), 


n a /s 


m«(T \WJpqIN? 

T? m 2 Tim p i 


[The quantity n a / & corresponds to the number of axions per comoving volume (oc n a R(t) 3 ). 
This is because so long as the expansion is adiabatic, the entropy density s scales as R(t)~ 3 . 
R(t) is the scale factor of the Universe.] 

The energy density today then is given by this constant ratio times the present entropy 
density (so — 7.04n 7 ~ 2970 cm” 3 ) times the axion mass. Remembering that T\ oc m 8 - 18 , 
we see that U a h 2 oc m” 1-18 , where the unusual power of the mass traces to the way in 
which the axion mass turns on (see Fig. 5). 

When this calculation is done very carefully (anharmonic effects taken into account, 
the motion of the axion field integrated precisely, etc.) the following expression results for 
the axion’s contribution to the present energy density 43 

f l a h 2 = 0.85 x 10 ±0,4 A7 o ° o -7 (m a /10” 5 eV)” 1,18 
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where A200 = Aqcd/200 MeV, and the 10 ±0,4 factor reflects the theoretical uncertainties. 
In deriving this formula, it has been assumed that there has been no significant entropy 
production since the epoch of axion production (T ~ Ti). If there has been significant 
entropy production, say the entropy per comoving volume (5 = i2(f) 3 s) increased by a 
factor 7, then is reduced by the same factor of 7. 37 

Moreover, it has been assumed that the initial misalignment angle of the axion field is 
just the rma average of a flat distribution of initial values from 0 to 7r, that is, = irjs/Z. 
Assuming that the Universe never underwent inflation this is the reasonable thing to do: At 
the time the axion field began to oscillate the presently observable Universe was comprised 
of about 10 30 or so causally-distinct volumes, each of which should have an independently 
choosen value for 0i ; clearly, the value of 9\ relevant to computing the average mass density 
of axions in the Universe today is 

On the other hand, if the Universe underwent inflation, either after or during PQ sym- 
metry breaking, then the entire observable Universe should be within a single inflationary 
region (or bubble, if you prefer) within which 61 takes on the same value. That value is 
equally likely to be in any interval between 0 and 7r, i.e., 61 is just as likely to be between 
0.1 and 0.2, as it is to be between 1.5 and 1.6. [Of course, it is fair to say that the a priori 
probability of 6\ being in the interval [1,1.5] is 100 times greater than it being in the 
interval [0.01,0.015].] If one averages over all inflationary patches in the entire Universe 
one can say that the rma value of 0\ should be / \/Z. However, that tells us nothing about 
the initial value of © in our neck of the woods. In order to determine 6 \ , we would have to 
measure the mass density of axions and the axion mass! Thus it is clear that in this case 
the energy density in axions is not precisely determined. Putting in the 6 \ dependence in 
the energy density of relic, coherent axions we have 

n a h 2 = 0.13 x 10 ±0 - 4 A 2 ~ 0 0 o 7 /(5 1 )^ 1 2 (m a /10- 5 eV)- 118 

where the function f{6) accounts for anharmonic effects: f{6) is monotonically increasing 
and f(6\ = 0) = 1.0. 

Note that the theoretical uncertainties inherent in fi 0 are large: from particle physics 
a factor of 10 ±0,4 , and from cosmology a factor of h 2 — all told, easily a factor of 10. For 
canonical values and 6\ — w/\/3, an axion mass of about 10 5 eV or so corresponds to 
closure density in axions. 

The final mechanism for axion production is even more intriguing: axion production 
through the decay of axionic strings, and it was first discussed by Davis. 44 In the case 
that the Universe never inflated, the initial value of 0 not only uniformily samples the 
interval of [0,x], but also has non-trivial topology. That is, the initial mapping of 6\ to 
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our 3- dimensional space cannot in general be smoothly deformed to a uniform value of 0i 
throughout space. The topological entities which exist are axionic strings. Let me be a 
little more specific about their formation and consequences. 

In most axion models PQ symmetry breaking is effected by a complex scalar field, 
denoted by B, which carries PQ charge. During PQ SSB B acquires a vacuum expectation 
value: < |<r | >= /pq. However, the argument of B is left undetermined — it is the axion 
degree of freedom. Consider the configuration around some axis where far from the axis 
< >— /pq exp(t^>), and <j> is the angle around the axis. This configuration cannot be 

smoothly changed into the configuration where the argument of < B > is constant and 
corresponds to an axionic string. In the core of the string (i.e., along the axis) < |<7| > 
must necessarily vanish; and owing to the vacuum energy associated with < B >= 0, 
together with the gradient energy associated with < B > (|<7| changes from 0 in the core 
to /pq far from the core) the string has an energy per unit length associated with it of 
order p ~ /pQl n (/pQd), where d is the characteristic distance between axionic strings. 

[The energy per length of a single, isolated string is logarithmically divergent, a well- 
known feature of global strings. That the spontaneous breakdown of PQ symmetry should 
result in topologically stable string configurations is expected since Hi (17(1)) = Zn- A 
spontaneously broken 17(1) symmetry is the simplest model for a string; for a review of 
strings and their properties, see ref. 45.] 

If the Universe did not inflate, then after PQ SSB, the Universe should be filled with a 
network of axionic strings. Much is known about the evolution of string networks. 46 Very 
rapidly a scaling solution is reached where the energy density in string resides primarily 
in a few infinite portions of string per horizon volume and is given by 

„ /* 

^string ~ 

Here t is the age of the Universe. It is referred to as the scaling solution because the ratio 
of energy density in string to the total density of the Universe is constant and equal to 
~ G/x. 

That the string energy density should evolve in such a way is somewhat surprising: On 
naive grounds, owing to the conformal stretching of the string network by the expansion of 
the Universe, one would expect the energy density of a string network to scale as R(t)~ 2 . 
In a radiation-dominated Universe this would imply that String « < -1 ; i/this did occur the 
energy density in string would rapidly grow relative to the radiation density, and string 
would soon come to dominate the energy density of the Universe. However, this does not 
occur because of dissipation: the cutting up of long pieces of strings into loops by string 
self-intersection and the dissipation of the oscillation energy of loops into gravitational 
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waves or other forms of radiation. In the present case, the dominant form of dissipation is 
the radiation of axions! 

In each Hubble time (H -1 ~ t) essentially the full energy density of axionic string 
is converted into axions, so that the change in the number density of axions per entropy 
density is given by 

%/>) ~ 

where during the radiation-dominated epoch the Hubble constant H ~ t~ l ~ T 2 /m p i and 
a? is average energy per axion radiated. The total number of axions produced per comoving 
volume is obtained by integrating the above equation: 

n a dT 

M / 2" 

s J Tl ujm^ 

where T\ is the temperature at which the axion mass becomes comparable to the expansion 
rate. After this, the axion mass becomes significant, and the string network becomes a 
network of domain walls bounded by strings which quickly decays. 

In order to calculate the axion production via this mechanism we must know u>(t), the 
average energy of an axion produced by string dissipation at time t. This is where the 
discussion heats up: Davis 44 argues that the axions produced have the longest wavelengths 
that they could be expected to have, of order the horizon, u>(<) ~ f -1 ; whereas Haraxi 
and Sikivie 47 argue that there is a 1/k spectrum of axion energies, which leads to u> ~ 
ln(/pQ<)f -1 — a difference of a factor of ln(fpqi) or about 100! The number of axions 
produced by string decay (per comoving volume) is 

n a /s ~ [1 or ln(.fr>Qfi)];p^- 

■Lirripi 


The form of this expression should be familiar: up to a factor of [1 or ln(/pQ<;i)]/0 2 it 
is identical to the expression for the number of axions produced by the initial misalignment 
of the © angle! Depending upon the spectrum of axions produced by the decay of axionic 
strings, these axions either contribute a comparable number of axions, or 100 times as 
many axions, as the misalignment mechanism does. 

More precisely, Davis 44 claims that the string-produced population of axions con- 
tributes a mass density 

n a h 2 


Harari and Sikivie 47 would claim that the contribution of string-produced axions is about 
a factor of 100 less. The difference is crucial: if Davis is correct axions provide closure 
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density for a mass of ~ 10 -3 eV, rather than ~ 10 -5 eV. If Harari and Sikivie are correct, 
then string-produce axions only slightly increase the axion density over that due to the 
misalignment mechanism. 

In any case we see that the energy density of axions produced by coherent processes 
increases with decreasing axion mass. Based upon our knowledge of the present age of 
the Universe, that it is greater than 10 Gyr, Clh 2 must be less than about 1. In the 
non-inflationary case, this restricts the axion mass to be 

m a > 10-« eV 

if Harari and Sikivie are correct; and 

m a > 10~ 3 eV 

if Davis is correct. Because of the peculiar scaling of fl a with the axion mass, the cosmic 
density provides a lower limit to the axion mass. 

In the inflationary case there are no string-produced axions (as the value of is 
uniform throughout the observable Universe). Moreover, the bound to m a based upon the 
present mass density of axions depends upon 6 \ : 

m a > 10- fl ft/Or/vW-'eV 
These bounds are shown in Fig. 2. 

We see that if Harari and Sikivie are correct, or if the Universe never underwent 
inflation, there is a substantial window between the SN 1987A bound and the axion mass 
density bound. On the other hand, if Davis is correct and the Universe never inflated, that 
window collapses to a single value for the axion mass: m 0 = 10 -3 eV. 

To summarize, the mass density of relic, thermal axions can probably never approach 
closure; however, the mass density can be significant enough to make an axion of mass 
2-5 eV detectable through its radiative decays. The mass density of relic axions produced 
by misalignment, or by the decay of axionic strings, can be very significant. However, 
there are still substantial uncertainties in the calculation of the relic axion abundance: if 
the Universe inflated, the value of 9\\ and if the Universe did not inflate, the spectrum of 
string-produced axions. In either case, the additional uncertainties in the values of Aqcd 
and the Hubble constant further complicate matters, making a precise determination of 
the value of the axion mass which leads to fl a — 1 impossible at present. 

One thing is certain, if axions do indeed contribute (l a = 1, they behave as cold dark 
matter because of their intrinsically small velocities. 48 And therefore, there is every reason 
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to believe that axions provide the dark matter known to exist in our halo, and therefore 
have a local density of 43 

pa = Phalo — 0.3 GeV cm -3 
n 0 ~ 3 x 10" 13 (10" s eV/m 0 ) cm" 3 

This estimate of the local axion density should be accurate to about a factor of 2. I should 
also mention that cold dark matter with inflation-produced adiabatic density perturbations 
continues to be a very attractive scenario of structure formation. 49 

[In axion models where N > 1, there tire actually N distinct, degenerate vacua: 50 
0 = 27rn (n = 0,1 ,..., N — 1). This means that in the absence of inflation different 
regions of space will wind up in different minima (i.e., different values of n), and will be 
separated by axionic domain walls of surface density cr ~ m T / T /pQ. Such domain walls are 
cosmologically very bad, leading to a wall-dominated Universe. There are number of ways 
of avoiding this catastrophe, including inflation. 51 Consideration of axion domain walls 
doesn’t really lead to any constraint to the axion mass: if they exist, axions are precluded, 
and so we will not address that issue here.] 

VI. Summary: Windows of Opportunity 

A priori the axion window spans an 18 order-of-magnitude mass window: 10~ 12 eV — 
10 fl eV. Laboratory experiments have probed masses greater than about 10 4 eV without 
success. This of course includes the originally favored value, m a ~ 200 keV, corresponding 
to a PQ SSB scale equal to that of the electroweak scale. Through a variety of very clever 
and interesting astrophysical and cosmological arguments, all but about 3 of the remaining 
16 orders-of-magnitude have been probed. The evolution of red giant stars precludes an 
axion mass in the intervals: 10" 2 eV— 10® eV (DFSZ) and 2/[{E/N— 1.95)/0.72] eV — 10 5 eV 
(hadronic, or DFSZ axion with /? ~ ir/2). The duration of the neutrino burst associated 
with SN 1987A precludes an axion mass in the interval: 10” 3 eV — 2eV (for either type 
axion). Radiative decays of relic, thermal axions preclude axion masses of 5-30 eV. 

Very low mass axions are precluded by the relic density of axions produced by coherent 
processes in the early Universe. In the case that the Universe underwent inflation, masses 
smaller than about 10"® eV[^i/(7r/'\/3)] 1 ' 7 axe precluded. Unfortunately, the value of the 
initial misalignment angle 6\ in our inflationary bubble is not known, and could with 
equal a priori probability lie anywhere in the interval [0, 7r], If the Universe never inflated, 
then the rm» value of 0\ is just 7 t/\/ 3. However, in this case there are two mechanisms 
for axion production: the misalignment of the 0 angle, and the decay of axionic strings. 
Production from the first mechanism is straightforward to compute (although there are 
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residual uncertainties), and an axion mass smellier than ~ 10 -6 eV is precluded. There 
is still spirited debate about the number of axions produced by the second mechanism: 
Davis 44 claims 100 times that of the misedignment mechanism; while Harari and Sikivie 47 
claim that the number produced is comparable to that of the misalignment mechanism. If 
Davis is correct and the Universe never inflated, then an axion mass less than about 10 -3 
eV is precluded. 

At present we are left with two windows of opportunity: (2 eV-5 eV) (for hadronic 
axions only) and (10 -8 eV — 10 -3 eV), corresponding to PQ SSB scales of few x 10 7 GeV 
and 10 10 — 10 13 GeV respectively. As noted previously, 10 -8 eV is a “soft boundary;” if 
the Universe inflated, then the precise value depends upon 8\‘ 7 , and if the Universe never 
inflated, and Davis is correct the 10 -8 — 10 -3 eV window shrinks to the single value of 
10- 3 eV. 

The 2-5 eV window for hadronic axions will be probed by at least two experiments: a 
telescopic search for photon line radiation from the decays of relic axions; 52 and an exper- 
iment designed to detect axions emitted by our own sun, which will be sensitive to axions 
of mass 0.1 — 5 eV. 53 The other window is even more intriguing because it encompasses the 
value of the axion mass for which axions provide closure density. Historically, that value 
has been taken to about 10~ 5 eV; as we have discussed, there are many caveats. If the 
Universe inflated, about all we can say is that the value which yields fi 0 = 1 is less than 
about 10 -5 eV and depends upon 8] 7 . If the Universe didn’t inflate, then the situation 
is equally or even more complicated. Depending upon the spectrum of string-produced 
axions, the axion mass resulting in fl 0 = 1 could be anywhere between 10 -5 eV and 10~ 3 
eV. Sikivie 54 has proposed a very clever idea for detecting these cosmic axions, based 
on the conversion of cosmic axions to monoenergetic photons of energy m 0 in a strong 
magnetic field. Already experiments 55 ’ 58 axe underway to search for relic axions of mass 
10~ 5 eV or so using his idea, and results have already been reported for axion masses of 
0.4 — 1.5 x 10 -5 eV. 55 However, the potentially interesting mass range up to 10~ 3 eV is 
not yet being probed, and stands as a real target of opportunity! 57 

Does the axion actually exist? I don’t know. However, I am very certain that if it 
does, it will be found in the heavens and not on earth! 

This review of the axion window is by no means complete; for more extensive reviews 
of the axion, see refs. 58. This work was supported in part by the DoE (at Chicago). 
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Table 1: Summary of axion masses excluded by astrophysical arguments based upon 
stellar evolution. Here £ = ( E/N — 1.95)/0.72. 


OBJECT 

DFSZ 

HADRONIC 

Sun 

leV-lOkeV 

20eV/<- lOkeV 

Red Giants 

10 -3 eV — 200 keV 

2eV/C- 200keV 

SN 1987A 

lO- 3 eV — 2 eV 

10 -3 eV — 2eV 

White Dwarfs 

3 x lO" 3 eV - few keV 

no limit 

Neutron Stars 

>10- 3 eV (?) 

> 10- 3 eV (?) 
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Figure Captions 


Fig. 1 . Axion couplings to ordinary matter (electrons, nucleons, and photons) and the dom- 
inant axion emission processes in stars. For the DFSZ axion the dominant emission 
processes in main sequence stars, red giant stars and white dwarfs are: the Compton 
like and bremsstrahlung processes, the rates for both of which are proportional to the 
axion-electron coupling squared. For the hadronic axion the dominant emission pro- 
cess in these objects is the PrimakofF process, owing to the fact that the tree level 
axion-electron coupling is highly suppressed. In neutron stars the dominant emission 
process for both types of axions is nucleon-nucleon, axion bremsstrahlung. 

Fig. 2 Summary of the laboratory, astrophysical, and cosmological constraints to the axion 
mass, and the two remaining axion windows. Note the constraint based upon the 
cosmic density of string-produced axions is still very uncertain, as it depends upon the 
spectrum of string-produced axions (an issue which is still in dispute). If the Universe 
underwent inflation after, or during PQ symmetry breaking, then there are no string- 
produced axions and the cosmological limit based upon coherent production due to 
the initial misalignment of the axion field depends upon the initial misalig nm ent, angle 
(to the 1.7 power). 

Fig. 3 The charactersitic length of the predicted neutrino burst in the KII and IMB H 2 O 
Cerenkov detectors as a function of axion mass for three different neutron star cooling 
models which include axion emission (from ref. 31). The quantity A<(90%) (in sec) 
is the time required for the expected number of neutrino events to achieve 90% of its 
asymptotic value. Note that for an axion mass greater than ~ 10 -3 eV the duration 
of the neutrino bursts becomes significantly shorter than those observed (~ 6 sec for 
IMB and ~ 12 sec for KII), thereby precluding such a value for the axion mass. 

Fig. 4 The axion luminosity from the nascent neutron star associated with SN 1987A (based 
upon a simple analytic model 32 ) as a function of axion mass. For m a < 0.02 eV axions 
simply freely stream out, and Q a oc m 2 ; for m 0 > 0.02 eV axions interact so strongly 
that they become trapped and are radiated from an axion sphere (like neutrinos). In 
this regime Q a oc T* oc For an axion mass in the interval [10~ 3 eV, 2eV], the 

axion luminosity (more precisely, cooling rate) is unacceptably large, precluding such 
a mass. 

Fig. 5 Summary of the relic axion contribution to Clh 2 from the 3 production processes: 
thermal, 39 misalignment, 37 and axionic string decay. 38 Note that if the Universe un- 
derwent inflation after, or during PQ symmetry breaking, there would be no string- 
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produced axions and the production due to misalignment would be proportional to the 
initial misalignment angle squared. Also note that the contribution of string-produced 
axions depends crucially upon the spectrum of axions from string decay, an issue which 
is still being debated, and which leads to an uncertainty of a factor of ~ 100 in axion 
production by this mechanism. 

Fig. 6 High resolution spectrum of the night sky at Kitt Peak 42 and the axion line expected 
from axion decays in our galactic halo. 41 Note that an axion of mass 4 eV would have 
produced a very prominent line — more prominent than is seen. 
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